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Abstract 

The PHOBOS experiment at the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory 
is studying interactions of heavy nuclei at the largest energies available in the laboratory. The high multiplicity 
of particles created in heavy ion collisions makes precise vertex reconstruction possible using information from a 
spectrometer and a specialized vertex detector with relatively small acceptances. For lower multiplicity events, 
a large acceptance, single layer multiplicity detector is used and special algorithms are developed to reconstruct 
the vertex, resulting in high efficiency at the expense of poorer resolution. The algorithms used in the PHOBOS 
experiment and their performance are presented. 
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1. Introduction 

The PHOBOS experiment is studying Au+Au, 
Cu+Cu, d+Au and p+p collisions at the nucleon- 
nucleon center of mass energies from ^/snn = 19.6 
GeV to y/sN~N = 200 GeV (and ^/JWn = 410 GeV 
for protons). For an overview of results see [1], for 
the description of the detector see [2,3]. General 
properties of these collisions are measured using a 
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Fig. 1. The central part of the PHOBOS detector [2]. Only 
the silicon sensors of the octagonal multiplicity detector, 
the vertex detector and the eight first planes of the two 
spectrometer arms arc shown. In the PHOBOS coordinates 
system mentioned in the text, the z axis coincides with 
the beam line, the y axis points upward and the x axis is 
horizontal in the direction of spectrometer arm P. 



large acceptance multiplicity detector covering al- 
most the full solid angle. About 50% of the pro- 
duced particles are registered in the central, octag- 
onal part of multiplicity detector (Fig. 1), whereas 
those emitted at smaller angles with respect to the 
beam line are detected in the multiplicity ring de- 
tectors placed a few meters from the nominal in- 
teraction point (not shown in Fig. 1). Detailed in- 
formation on about 1% of the particles is obtained 
from a two-arm magnetic spectrometer. The posi- 
tion of the primary interaction is precisely recon- 
structed using a specialized vertex detector, which 
registers about 5% of the primary particles. At 
RHIC, the beams of ions collide at zero degree, 
thus the primary vertices are distributed in a wide 
z range along the beam line (up to ±2 m from the 
nominal interaction point) and only about ±0.1 cm 
in the perpendicular direction. The relatively small 
acceptance of the vertex detector is sufficient for 
vertex reconstruction for events with large multi- 
plicities, but does not allow to extend the anal- 
ysis to more peripheral Au+Au collisions or to 
the interactions of lighter systems, especially p+p. 
Therefore, in order to improve the vertex determi- 
nation efficiency and accuracy, several algorithms, 
presented in this paper, are needed using comple- 
mentary information from different detector sub- 
systems. 



2. Vertex reconstruction algorithms 

The scintillator counters of the trigger system 
(not shown in Fig. 1) deliver the first information 
on the position of the primary vertex. Two larger 
detectors, paddle counters, are placed symmetri- 
cally at z = ±3.2 m from the nominal interaction 
point, two smaller TO counters are mounted around 
the beam pipe at z = ±5.3 m for most measure- 
ments (in p+p collisions they are moved closer, 
to z = -1.3 m and z = +2.5 m). The coincidence 
of the signals from these detector pairs on differ- 
ent z-sides is used for triggering, but as the arrival 
times of produced particles are measured, the ver- 
tex z position can be derived from the registered 
time difference. The error in the vertex position 
determination (about 5 cm) is small enough for an 
on-line vertex range selection, but it is not suffi- 
cient for data analyzes. 

The vertex position in 3 dimensions (X v , Y v , Z v ) 
can only be obtained from tracks reconstructed in 
the PHOBOS spectrometer. It consists of 16 lay- 
ers of silicon pad and strip sensors places inside 
the conventional 2 Tesla magnet. The design of the 
magnet ensures that the first 8 layers of spectrom- 
eter silicon sensors are in an area where the field is 
very weak. This allows easier reconstruction of the 
first segment of the particle trajectory as a straight 
line, which is then normally joined with a curved 
part of the track found in the remaining layers. 
The acceptance for the straight line segments is 
larger than that for the fully reconstructed tracks 
and the information on the curved part does not 
improve track extrapolation to the vertex, thus the 
straight lines are used in two algorithms of vertex 
reconstruction. 

In the spectrometer 3D method, the vertex po- 
sition is defined as the point for which the sum 
of distances from the extrapolated tracks is mini- 
mal. In an iterative procedure, tracks too distant 
from the calculated vertex (eg. secondary parti- 
cles) are removed and the vertex is recalculated. 
The spectrometer 2D+1D algorithm enters the 
points of closest approach of all track pairs into 
2-dimcnsional (x — z) and 1-dimcnsional (y) his- 
tograms. The coordinates of the maxima in both 
histograms define the approximate vertex posi- 
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tion. The final vertex position is determined by 
the mean in x, y and z of only those points com- 
patible with the beam orbit and close enough to 
the approximate vertex. The sizes of the pads in 
the silicon sensors used for vertex reconstruction 
are relatively large (lxl mm 2 and 0.427x6 mm 2 ). 
This, together with small acceptance of the spec- 
trometer, is responsible for the limited precision 
of the spectrometer methods. 




Fig. 2. Determination of the vertex position using all com- 
bination of hits registered in the vertex detector. Dashed 
lines represent false combinations, continuous lines primary 
particles. The histogram of the number of tracks extrapo- 
lated to a fixed y value (at the bottom of the picture) is 
used to find the position of the vertex. 

The PHOBOS vertex detector consists of two 
parts placed above and below the beam pipe, 25 cm 
long in z, each with 2 layers of silicon sensors at 
about 5.7 cm and 11.8 cm from the beam. The sen- 
sors have strips 473 /im wide and 1.2 or 2.4 cm long, 
perpendicular to the beam direction. The strips in 
the more distant outer layers are two times longer 
to cover the same azimuthal angle range. The de- 
sign of the vertex detector limits the vertex re- 
construction to y and z positions. In the vertex 
search all combinations of hits in the two layers are 



treated as track candidates and their extrapolated 
z position at a fixed y is histogrammed (Fig. 2). 
Z v is defined by the position of the maximum of 
such a histogram. This procedure is repeated for 
several values of fixed y. The value of y for which 
the maximum of the histogram is the highest is 
used as final Y v . To accept the maximum at least 
three track candidates contributing to it are re- 
quired, since without this cut coincidences of two 
false track candidates are selected in low multiplic- 
ity events more frequently than the true vertex. 
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Fig. 3. The correlation of the difference between the po- 
sition of the hit and the position of the primary vertex, 
Zhit ~ Z v , and the registered energy loss AE norm in the 
PHOBOS octagonal multiplicity detector. To correct for 
the different sensor thicknesses, the measured energy loss 
is divided by the mean energy of a minimum ionizing par- 
ticle traversing the sensor perpendicularly. 

In PHOBOS, the largest fraction of primary 
particles is registered in the octagonal multiplic- 
ity detector consisting of 92 silicon sensors with 
0.27x0.86 cm 2 pads. Their thickness varies from 
0.305 mm to 0.338 mm. The sensors are placed 
around the beam pipe at radii of about 4.5 cm 
and cover the z range from -55 cm to +55 cm. In 
this detector primary particles traverse only one 
silicon sensor and do not leave a track directly 
pointing to the vertex. However, the length of 
their trajectory inside the silicon (and thus the 
deposited ionization energy, AE) depends on the 
angle at which they enter the sensor; for primary 
particles approximately the emission angle. In this 
way the measured energy loss is correlated with 
the distance of the hit from the vertex (Fig. 3). 
Most of the hits belong to one of the two ridges, 
symmetric with respect to Zhu — Z v = 0. The hits 
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left by particles with very low momenta, which 
deposit much larger energy than the minimum 
ionizing particles, are visible at Zhit — Z v close to 
zero as a tail of large AE values. Because of the 
positive/negative Zhit — Z v ambiguity, the value 
of AE from a single hit can be used only to esti- 
mate the distance to the vertex or two alternative 
vertex positions. To determine the z vertex po- 
sition slices with fixed AE values are calculated 
from the 2-dimensional histogram shown in Fig. 3. 
Several examples of these histograms are shown 
in Fig. 4. Using such histograms the probability 
density function, P(AE, \Zhit — Z v \), for a particle 
that loses the energy AE at the position Zhit to 
originate from the vertex at Z v is calculated. This 
function is non-zero in two z ranges, symmetric 
with respect to Zhit- For most registered primary 
particles the width of a single range is 3-6 cm, and 
useful information is provided by the hits at the 
distance \Zhit — Z v \ up to 100 cm. For a single hit 
there is no way to select the correct range, and 
one needs an overlap of several ranges from more 
hits to find the correct vertex position. The z po- 
sition of the vertex is determined by maximizing 
the product of the probabilities of all hits in the 
octagonal multiplicity detector. 




Fig. 4. Distributions of the distances between the primary 
vertex and the hits for several different values of the nor- 
malized energy loss AE norm . The histograms (from left to 
right) arc obtained for AE norm equal 1.2, 3, 6, 15 and 30. 




1- 
O.S r 
0.6- 
0.4 ^ 
0.2- 



b) 




-40 



-20 



20 



40 



60 




-%0 



-40 



-20 



20 



40 



60 



Fig. 5. The vertex reconstruction efficiency for Au+Au 
collisions at y/s NN = 200 GeV defined as fraction of events 
for which the distance between the reconstructed and true 
vertices does not exceed 0.5 cm, 0.1 cm and 3 cm for the 
three methods: spectrometer 2D+1D (a), vertex detector 
(b) and octagon (c). The additional events with incorrectly 
reconstructed vertex arc represented by the shaded area. 
The maximal allowed distances between the reconstructed 
and true vertices are adjusted to the resolution of the 
method. 

3. Vertex reconstruction quality 

In the analysis of the vertex reconstruction qual- 
ity we will focus first on Au+Au collisions at the 
highest RHIC energy, y^s^ — 200 GeV, using 
Monte Carlo simulations of the PHOBOS detector 
response based on GEANT 3.21 package [4]. 

The PHOBOS vertex detector is optimized for 
measurements of collisions occurring close to the 
nominal interaction point (Z V =Q), but the ac- 
ceptance of some of the reconstruction methods 
extends over much larger Z v range (Fig. 5). The 
octagon method is the most efficient (about 95% 
for the inclusive sample of Au+Au collisions) and 
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for all events, in which the vertex can also be re- 
constructed using other, less efficient methods, it 
returns the correct position. All methods are op- 
timized to reconstruct the vertex even for events 
with small multiplicities, thus for a fraction of 
events they fail to find the correct vertex position. 
In the case of the spectrometer 2D+1D method 
that results in a larger than typical error, while 
the vertex detector method quite frequently finds 
a vertex close to z=Q in the events, where the 
true Z v is outside the acceptance of this detector. 
Fortunately, for the later case the octagon method 
gives a different vertex position and, thus, such 
events can be eliminated from the analysis. 
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Table 1 

Vertex position accuracy ([cm]) and efficiency for Au+Au 
collisions at ^/s NN = 200 GeV in the vertex range lim- 
ited to \Z V \ <10 cm. To the sample of the most central 
events belong these with the largest values of the number 
of charged primary particles hitting the octagonal multi- 
plicity detector, to the peripheral sample these with the 
smallest multiplicities. For the most central events the effi- 
ciency of all methods is 100% (not shown). In order to re- 
duce the contribution of the tails of large errors the quoted 
resolution is taken from the fit of a Gaussian rather than 
a simple RMS of the ^ me t?j,od — %v distribution. 

The reconstruction error depends strongly on 
the multiplicity of the events, which becomes ev- 
ident if from the sample of Au+Au collisions the 
most central and the most peripheral events are 
selected and separately analyzed (Table 1). In the 
first subsample all methods are 100% efficient and 
the vertex detector method is the most accurate in 
determining Y v and Z v . For peripheral events only 
the octagon method has a large efficiency (85%), 
whereas both spectrometer methods are inefficient 
(4-7%). The error of the vertex position increases 
for all methods, about a factor of 2 for the octagon 
and more than a factor of 10 for the spectrometer 
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Fig. 6. The vertex reconstruction efficiency as a function 
of the number of primary charged particles present in the 
subdetector used by spectrometer 2D+1D, vertex detector 
and octagon methods. 

3D method. The vertex detector method is able to 
reconstruct many more events than the spectrom- 
eter methods and still gives the most accurate ver- 
tex position. 

In order to compare reconstruction results from 
Monte Carlo simulations and real data, we also cal- 
culate the values of a{Z met hod - Z verte x -method) for 
spectrometer 3D, spectrometer 2D '+ ID and octagon 
methods using events, in which all methods have 
found the vertex (about 50% of all Au+Au events). 
In the case of the spectrometer methods the errors 
agree very well. The error for the octagon method is 
about 20% smaller in the real data than in the sim- 
ulations. This can be explained by an inaccurate 
description of the energy loss distribution or by a 
difference in the contribution of background par- 
ticles in the GEANT simulations. Especially the 
slightly longer tails of the energy loss distributions 
in the simulations can lead to the observed larger 
error in the calculations of the octagon vertex po- 
sition. 

The efficiencies of the vertex reconstruction 
methods used in the PHOBOS experiment reflect 
the acceptance of the subdetectors used by the 
reconstruction algorithm, mentioned in the Intro- 
duction. To compare the "reconstruction power" 
of the methods using one, two or several layers of 
silicon we calculate their efficiencies as the func- 
tion of the number of primary particles which 
point to the vertex and are registered in the appro- 
priate subdetector (Fig. 6). The high quality of the 
information from reconstructed tracks causes that 
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Table 2 

The efficiency of the vertex reconstruction algorithms for 
different collision types. The vertex range was limited to 
\Z V \ <10 cm. The mean number of charged primary par- 
ticles with hits in any subdetector is given as an approxi- 
mation of the number of produced charged particles in the 
analyzed samples. 



spectrometer 2D+1D method can find the vertex 
in more than 60% of events, in which only two 
primary particles traverse the spectrometer, and 
from 4 such tracks it is over 90% efficient. Inter- 
estingly, sometimes the vertex is also found in the 
events with only one primary particle, but then at 
least one non-primary particle, for example from 
a decay very close to the vertex, helped in the 
reconstruction. For the vertex detector method, 
as enforced by the reconstruction requirements, 
the minimal number of primary particles is 3 and 
100% efficiency is reached for 15 such tracks. The 
octagon method needs at least 5 primary particles 
and is more than 80% efficient if 10 primary par- 
ticles are registered in the octagonal multiplicity 
detector. 

Finally, it is worth to mention the overall effi- 
ciency of the vertex reconstruction methods for dif- 
ferent types of interactions measured by the PHO- 
BOS experiment (Table 2). For collisions of heav- 
ier nuclei (Au+Au or Cu+Cu) the efficiencies of 
the precise spectrometer and vertex detector meth- 
ods are close to 50% or higher. In the case that 
deuteron is one of the projectiles, the efficiency 
drops faster for the spectrometer method, and be- 
comes very small for both methods in the case of 
p+p collisions. The octagon method has an effi- 
ciency of more than 86% for nucleus+nucleus and 
more than 50% for p+p interactions. 



4. Summary 

The PHOBOS experiment utilizes several vertex 
reconstruction algorithms, which use complemen- 
tary information from different subdetectors. Ac- 
cording to the design goals, a very precise vertex 
resolution is achieved in high multiplicity heavy 
ion collisions, but the limited acceptance of the 
spectrometer and vertex detector leads to low effi- 
ciency in peripheral Au+Au collisions, d+Au and 
p+p interactions. To overcome this problem a new, 
very efficient method of vertex reconstruction us- 
ing hits from a single layer octagonal multiplicity 
detector is introduced. Even for low multiplicity 
events it efficiently finds the vertex with an error 
of about 1 cm. This method can either be used to 
verify doubtful results of the other algorithms or 
directly in studies, for which such precision is suffi- 
cient. Such procedure for testing the vertices from 
different methods is implemented in the PHOBOS 
experiment software. It rejects wrong reconstruc- 
tion results and calculates the most accurate ver- 
tex position. 
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